The growth of linear perturbations in generic defect models for structure formation 
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We study the growth of linear perturbations induced by a generic causal scaling source as a function 
of the cosmological parameters h, f2™ and 0.q . We show that for wavenumbers k > 0.01 h Mpc -1 the 
spectrum of density and velocity perturbations scale in a similar way to that found in inflationary 
models with primordial perturbations. We show that this result is independent of the more or less 
incoherent nature of the source, the small scale power spectrum of the source and of deviations from 
scaling which naturally occur at late times if Q$ 7^ 1. 
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A. Introduction 

At present there are two main competing paradigms 
for the origin of the large scale structure of the Universe: 
topological defects and inflation. In defect scenarios 
density perturbations are causally seeded by defect evo- 
lution on increasingly larger scales while in inflationary 
models the causal structure of the universe is altered and 
correlations can be established on scales much larger than 
the Hubble volume. Inflationary perturbations are pro- 
duced at a remote epoch and their evolution is linear 
until late times. They are passive and coherent. In con- 
trast the evolution of topological defects is highly non- 
linear and can only be fully described resorting to large 
numerical simulations |^,|3| . Defect perturbations are ac- 
tive and more or less incoherent depending on the source. 
Despite the differences between them both scenarios pro- 
duce a nearly scale invariant spectrum of density pertur- 
bations. Hybrid models featuring both inflationary and 
defect perturbations have also been proposed [Q ||. 

Recently there has been substantial progress in the un- 
derstanding of defect models as seeds for the large scale 
structure and cosmic microwave background anisotropics 
|0-|f2|. The matter and microwave anisotropy power 
spectrum induced by local and global strings, global 
monopoles and global textures have been calculated 
and compared with current observations. Although 
the results seem to indicate some disagreement between 
standard defect models and observations (in particu- 
lar in what concerns the cosmic microwave background 
anisotropies), some issues still remain to be investigated 
in more detail. This is particularly true for local strings 
in which the exact contribution of cosmic string loops 
and gravitational radiation is still uncertain Jf^Jf^] ■ 

Given that current observational data strongly suggest 
that the density of the Universe is sub-critical (favoring 
a non-zero cosmological constant) [[l5 16 it is very im- 
portant that structure formation analysis performed for 
defect models can be generalized for any realization of the 
cosmological parameters h, and fi A . This has been 
done in particular for cosmic strings 0JT^,|I^| but the va- 
lidity of the approximations used have not been investi- 



gated for generic defect structure formation scenarios. In 
this letter we describe the dependence of the power spec- 
trum of density and velocity perturbations, induced by a 
generic causal scaling source, on the cosmological param- 
eters h, fig 3 and fi A . We show that for k > 0.01 /iMpc -1 
these scale in a very similar way to that found in inflation- 
ary models with primordial fluctuations. We also study 
the deviations from scaling which naturally occur at late 
times if fi m starts deviating from unity. This is modeled 
by switching off the source when f2 m reaches a certain 
fixed critical value f2™ < 1. We demonstrate that our 
results hold independently of the more or less incoherent 
nature of the source, the small scale power spectrum of 
the source and of deviations from scaling which naturally 
occur at late times if i7 n 1 . Finally, we generalize 
our results for an arbitrary red-shift. 



B. The growth of linear perturbations 

The dynamical equation which describes the evolu- 
tion of the scale factor a in a Friedmann-Robertson- 
Walker (FRW) universe containing both cold dark matter 
(CDM) and radiation fluids can be written as: 
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where TL = a/ a (a dot represents a derivative with re- 
spect to r = ri/ri e q and r) is the conformal time), fl m , 
n r and fi A express the densities in CDM, radiation and 
cosmological constant as fractions of the critical density 
and Q k = 1 - fi m - Q 1 ' - fi A . The subscript '0' means 
that the quantities are to be evaluated today and the 
subscript 'eq' means that the quantities are to be evalu- 
ated at the time of equal matter and radiation. We have 
also assumed that ao = 1 and H(tq) = r] cq . 

The evolution of radiation and CDM density fluctua- 
tions (<5 r = Sp r /p r and S m = 6p m /p m respectively) in the 
synchronous gauge is given by ]20| : 

<5k + ns* - \n 2 {n m s^ + 2n r s£) = x>(k, r), (2) 



1 



si 



— k^Slr 



4 •• 

3 5k 



o, 



(3) 



where Si 



^k" 1 ^ 2 ' ^k' m 1S ^°- e Fourier transform of 



<5 r,m , k = |k| is given in units of (c?7cq) -1 , 2?(k, r) = 
47rfc 3 / 2 G8+, 6+ = 9oo + 6 M and 9 Q/ g is the energy- 
momentum tensor of the external source. 

At early times when a < min((r2™/r2o )i, Slg 1 /S1q) tnc 
scale factor is a quadratic function of the conformal time: 



o(t) = a oq (2(V2 - l)r + (3 - 2V2)t 2 ) 



and 
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with a cq cx (f2™/i 2 ) _1 . For a scaling source no other 
length scale apart from the horizon scale is involved and 
so 



with 



P(k,r)=0 2 fc 3 / 2 r- 1 /^ (k;T ) j 



|P(k,r)| 2 ) =^fc 3 r- 1 ^ av (fcr), 
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due to homogeneity and isotropy. Here 4>o is the defect 
symmetry breaking scale. 

We can see from equations that for a scal- 

ing source the evolution of density perturbations on a 
wavenumber k up to a conformal time r/ satisfying 
a{r f ) < min((fig7n£)i,ng7fi k ) does not depend on the 
cosmological parameters h, iljf and Qq (note that k is 
given in units of (c77 cq ) _1 cx fl^h 2 Mpc" 1 and r = n/r] eq ). 
We note that this conclusion still holds in the more re- 
alistic case in which the right hand side of equation (||) 
is multiplied by a function of r to account for the de- 
viation from scaling which naturally occurs during the 
radiation-matter transition. Consequently, we conclude 
that in order to study the dependence of the normaliza- 
tion of the power spectrum of density fluctuations, in- 
duced by a scaling source, on the cosmological parame- 
ters h, f2™ and we only need to study the evolution 
of perturbations from a conformal time Ti onwards (with 
ftjj/fig 1 < a{n) < min((^ n /^)i,^ n /^)). i n this 
case equation (||) reduces to: 



$k + HS£ -rLt&SP/a = Z>(k, r) 



(8) 



We solve this equation numerically assuming a simple 
form for the scaling source: 

T(k, r) = T{kr) = C^9(/3fcr - 1)6(1 - 0*kr), (9) 

with > 1, C = constant and 9 (x) is the step func- 
tion (9(x) = if x < and Q(x) = 1 if x > 0). The 
initial conditions in equation (||) are S^(n) = and 



Sk( T i) = f° r fi^Ti < 1. This source incorporates the 
most important features of causal scaling models relevant 
to our study. The parameter (3 is related to the compen- 
sation scale while /3» is a small scale cut-off which models 
the deviation of the power spectrum of the source from a 
white-noise spectrum on small scales. For the moment we 
are assuming that the source is coherent but in the next 
section we will see that this does not affect our results. 

Different values of £7™ and £!q lead to different lin- 
ear growth factors from early times to the present. For 
primordial perturbations the quantity £l™h 2 g(Q™ , Hq ) 
where 



7(«8\fi£) = 



5f^/2 



(ngov* - nft + (i + ngV2)(i + n£/ 7 o) 



(10) 



provides a very good fit to the growth factor of den- 
sity perturbations from early times (SIq/Oq 1 <C a(ri) <C 
rrun((Q^/fi£)i,ft|?/Q k )) to the present Q ( see also 
p2|). This fit was shown to be good to a few percent 
for 0.1 < O™ < 1 and < fig < 1 and is normalized to 
unity for Og 1 = 1 and f$ = 0. 

A quantity directly related to the amplitude of mat- 
ter perturbations is their rate of growth which may be 
described by : 



f = 



dlnS 
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For primordial perturbations a very accurate approxima- 
tion to the present rate of growth of matter perturbations 
is given by p§ 



/(fi- fi£) = (fi-)°' f ' 



70 



(12) 



In the next section we will show how these results can be 
generalized for generic defect models of structure forma- 
tion. 



C. Results and discussion 

The power spectrum induced by a scaling source for 
arbitrary values of the cosmological parameters h, f2™ 
and Oq can be written as 

s k (h, og\ n° A ) = a ■ 5 2 (n™ n£) . ^h 2 ) 2 ■ s k (i, 1, o), 



(13) 



where 



s k (h,n?,s#) = (\s k (h,nz,n*,T )\ A ), (w) 

A = <$(Off,n£)/($(l,0) and k is given in units of 
(c7? eq ) _1 (or equivalently, in units of f^/i 2 Mpc -1 ). The 
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FIG. 1. Plot of e g and £/ as a function of Sl m for open 
universes (solid lines) and flat universes (dashed-lines) in the 
case with (3/(3, -> oo for k = 0.05, 0.01 and 0.002 /iMpc" 1 in 
ascending order. 
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FIG. 2. Plot of e g and e/ as a function of fi m for an open 
universe in the case with (3/f3 t — > oo for k — 0.01 h Mpc -1 
and SI™ = 0,0.5,0.7 (solid, dot-dashed and dotted lines, re- 
spectively) . 



factor A(n$,tt$) can be normalized using, for example, 
the cosmic microwave background observations. 
It is also useful to define 



rfln|S k | 
din a 



fk 



in a similar way to equation ([LI 
The functions 



and 



e g (Je) = \l-g k /g\ x 100, 



e f (k) = \1 - h/f\ x 100, 



(15) 



(16) 



(17) 



measure the deviation of fk and gk from / and g, given 
respectively by equations ( |l0| ) and (|lj), evaluated at the 
present time. 

In figure [l] we plot e g and e/ as a function of Sl m for 
open universes (solid lines) and flat universes (dashed- 
lines) in the case with (3/ (3„ — > oo. We can see that gk = g 
and fk = f for k > 0.01 /iMpc -1 to high accuracy. We 
note that the case with /3//3* — > oo is the worst possible 
scenario as far as the smallness of e g and ef is concerned. 
For < 2, for example, gt ~ g and /fe ~ / to better 

than 5% for fc > 0.003 /iMpc" 1 . 

Here we have chosen a specific value for /3 (/? = 1/ vl8) 
valid for the matter era and motivated by the results of 
ref. p4|. The change in (3 which occurs at late times if 
Sl^ 7^ 1, does not alter our results. We have also verified 



that our results are not modified if we replace our sharp 
cut-off on large scales in equation (|^) by the more realistic 
cut-off 1/(1 + (fc/3r) 2 ) 2 . 

We also investigate the effect of deviations from scal- 
ing which may occur at late times by switching off the 
source for Q m < O™ < 1 . In figure ^ we plot e g and e / 
as a function of SI™ for open universes in the case with 
(3/13., oo for k = 0.01 /iMpc" 1 and Sl m = 0,0.5,0.7. 
We can see that our results are not greatly influenced by 
the late time behavior of the source. This is specially true 
if /?//?* < 5 or for flat universes, in which case e g j < 0.05 
for S1J 1 < 0.9 and fc > 0.01 h Mpc" 1 . 

Using the linearized continuity equation 



-«k • v, 



(18) 



where v is the peculiar velocity, it can also be inferred 
from our results that the velocity spectrum at the present 
time can also be rescaled for arbitrary values of the cos- 
mological parameters h, SI™ and S1q in a simple and ac- 
curate way: 

v k (Sl m , S$) = h- /(Sl m , n£) • jW, ) ■ v k (l, 0), (19) 

for k > 0.01 /iMpc -1 . 

We have also verified that our results hold for an inco- 
herent source satisfying : 



(e + (k,r)e + (-k,r + Ar)> =0, 



(20) 



for At > t c where r c is assumed to be smaller than one 
Hubble time and 



3 



(e + (k,r)e + (-k,r + T'))dr\ (21) 



is given by equation (Q). 

Finnaly we note that our results may be generalized 
for an arbitrary red-shift z = 1 /a — 1 : 

q fh m a s s k (h,n$,n£) g 2 (» m ,» A ) r99 , 

bk\h,il ,11 ,Z) = — ^ _9^m nA \ 1 <. 22 ) 



for fcr > 0.02. 

Although we have assumed the dark matter to be cold 
our results are also valid for hot dark matter due to scal- 
ing of the neutrino free-streaming length with fl^h 2 (see 
for example p5fl). 



D. Conclusion 

In this letter we have shown that the spectrum of den- 
sity and velocity perturbations induced by a causal scal- 
ing source for a particular realization of the cosmolog- 
ical parameters h, f2™ and A can be generalized in a 
very simple and accurate way for any reasonable choice 
of these parameters on scales k > 0.01 /iMpc^ 1 . This 
rescaling greatly simplifies the study of structure forma- 
tion with topological defects because it allows for an easy 
generalization, on all scales of cosmological interest, of 
the results of structure formation simulations induced by 
causal scaling sources for any reasonable values of the 
comological parameters h, 0™ and fi A . We have verified 
that our results are generic in the sense that they hold 
independently of the more or less incoherent nature of 
the source, the small scale power spectrum of the source 
and of deviations from scaling which may occur at late 
times. 
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